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Abstract 

In the present work a closed system of kinetic equations is obtained for the de- 
scription of the vacuum creation of an electron - positron plasma and secondary 
photons due to a strong laser field. An estimate for the photon energy distribu- 
tion is obtained. In the Markovian approximation the photon distribution has 
a 1/k spectrum (flicker noise) . 



Introduction 

Up to now, the Schwinger effect pQ has resisted an experimental verification. 
This is basically due to the huge critical field strength E c which could not yet be 
reached in the laboratory. Recently, main attention was devoted to theoretical 
studies of pair creation by time- varying electric fields [U El H] where sufficiently 
strong electric fields can be achieved at modern high-intensity laser facilities. It 
has been shown [2j[3j[4l[5] that pair creation by a single laser pulse with E <C E c 
could hardly be observed. More optimistic results have been obtained for a X- 
ray free electron lasers [7J [H] and for counter-propagating beams of optical 
lasers [HI HH1 EH] - It is obvious, that for subcritical fields E -C E c the electron 
- positron excitations have quasiparticle character and S - matrix methods can 
not be applied [12] and existing estimates [8] are not reliable. An adequate 
method is the kinetic theory. 

We will construct here the system of kinetic equations for a selfconsistent 
description of the electron - positron - photon system generated from the vacuum 
by a time-dependent electric field. As a first step, we will consider the one- 
photon annihilation process only, which in the presence of a strong field is not 
forbidden [15] . 

Below we will assume an external electric field A^(t) — (0, A oxt (i)) that is 
spatially homogeneous. In the case of a strong external electric field A cx t(£) also 
some internal field A; nt (t) will be generated. The total acting field will be equal 
to A(t) = Aj n t(t) + Aext(t) and this field is quasiclassical. Some fluctuations 
of the electromagnetic field can arise against this background. They can be 
interpreted as photon excitations. These photons, in principle, can be registered 
far from the active zone. 



1 E-mail: smol@sgu.ru 



1 Electron - positron sector 

In general, the complete system of equations for a selfconsistent description of 
the electron - positron - photon plasma consist of: 1) the KEs for the electron 
and positron quasiparticle components with distribution functions f e . P (p,t) in 
the presence of a total electric field E = — A(t), 2) the KE for the photon com- 
ponent, and 3) the Maxwell equation for the internal field Ai nt (t). We assume 
the electroneutrality condition f e (p,t) = f c {—p,t) = f(p,t) to be fulfilled. 

We start from the standard QED Lagrangian£ = £q+C taking into account 
in Cq the interaction with a quasiclassical (A^t)) and in £ with a quantized 
(Afj,(t)) electromagnetic field, 

C = ^{^D^-(D^)^}-m^ , £' = -e^l M tf, (1) 

where = + ieA^(t). It is assumed that the intensity of the quantized 
field is rather weak so that there is no backreaction influence on the state of 
the system. In other words, the electron - positron system plays the role of a 
photon source only. The in - vacuum \in) — |) is defined such that (A^x)) = 0. 
Below we will not consider the backreaction problem because for subcritical 
fields E <C E c the internal field is negligible. 

The kinetics of electron - positron vacuum pair creation due to a linearly 
polarized electric field was studied in a large number of works. See, e.g., Refs. 
[141 115] and the works quoted therein. The corresponding generalization to the 
case of an arbitrarily polarized time dependent electric field was obtained in 
[TBI riTl [T8] . The oscillator representation in its different realizations [IH [TS] 
leads to the nonstationary orthonormalized spinor basis 

uf(p,t)=A(p) [w + ,0,p 3 ,p_] , ut(p,t)=A(p) [0,w + ,p + ,-p 3 ] , 
v+(-p,t) = A(p) [- P 3 ,- P -,u+,0] , v+(-p,t) = A(p) [- P+ , P 3 ,0,uj + ] , (2) 

where ui(p, t) = V m 2 + P 2 = u, P — p — eA, p± — p 1 ± ip 2 , lo + = uj + m and 
A(p) = [2WW+]- 1 / 2 . 

Finally, the Dirac equation in the presence of an external quasiparticle field 
Aext(i) generates the Heisenberg-like equations of motion for the construction 
operators 

a(p,t) = -U (1) (p,t)a(p,t) - U (2 )(p,t)b + (-p,t) - icv(p,t)a(p,t) , 
b{-p,t)=b(-p,t)U {1) (p,t)+a + (p,t)U {2) {p,t)-iuj(p,t)b(-p,t) , (3) 

where in the representation ([2]) the matrices are U^(p, t) = iuja{pE,]cr = iU^a^ 
and U( 2 )(p,t) = q[P(PE) - Ecjoj + ]<t with a = e/ (2lo 2 uj + ). 

KEs for the electron - positron component of the plasma follow from the 
equations of motion ([3]) and the definitions of the electron and positron distri- 
bution functions in the instantaneous representation 

f a0 (p,t) = {a+(p,t)a a {p,t)) , /^( P ,t) = (bp(-p,t)b+(-p,t)) (4) 

and also the two additional functions 

/$(p,i) = (a+(p,t)b+(-p,t)) , f l J(p,t) = (b (-p,t)a a (p,t)) (5) 
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describing vacuum polarization. The system of KEs is then [T51 [T71 [TH] 
f=[f,U(i)}-{U (2) f i+) +f { - ) U {2) ), 

f c =[f c ,U (1) } + (f^U {2) +U (2) f^), (6) 
/(+) = [/(+), + (u {2) f - f c U (2) ) + 2*w/M 
/(-) = [/(-), U(1) ] + (/C/ (2) - U {2) f c ) - 2tu/(-). 

If the standard decomposition in the basis of Pauli matrices is used f = fo + 
fk&k, where f = Tr{/}/2, and = Tr{/<7fc}/2, the system KEs © can be 
rewritten in the spin representation. Some of the simplest applications of the 
corresponding system of KEs can be found in [T^l [TTl [TS] . 

2 Photon sector 

In order to construct the photon kinetics it is necessary to derive the correspond- 
ing generalization of the quasiparticle formalism developed in Sect. 1, using it 
as a nonperturbative basis. The interaction with the quantized electromagnetic 
field A^(x) in the fermion sector of the theory is introduced by means of the 
substitution Ho — ► Ho + H' in the Heisenberg-like equation of motion Q , 

d(p,i) + U (1) (p,t)a(p,t) + U {2) (p,t)b+(-p,t) = -i[a{p,t),H +H'}, (7) 

where Hq is the Hamiltonian of the fermion field in the quasiparticle represen- 
tation and H' is the usual Hamiltonian of interaction with the quantized field. 
The time dependence of H'(t) makes manifest the nonstationarity of the system 
that is also reflected in the decomposition of the field operators ^ + in the 
nonstationary basis (f5]). The same source (external field) induces the nonsta- 
tionarity of the quantized electromagnetic field. However, this does not alter 
the mass shell of the photon field, k 2 = (in contrast to electron - positron 
field, w(p,i)), so that the standard decomposition is valid, 

A,{x) = (2^)-V 2 | ^LA,(k,t)e-^, (8) 

where (k, t) = A ( + ] (k, t) + A<T> (-k, t). 

The system of Heisenberg-like equations of motion taking into account the 
photon subsystem can be written now in explicit form. For example, 

iA^^t) - TfcA (±) (k,i) T e(2n)- 3 /^ f d 3 Pl d 3 p 2 S( Pl - p 2 T k) 

|a+ (pi , t) [uu] r (pi , p 2 , t)a(p 2 ,t)+a + (p x , t) [uv] r (pi , p 2 , t)b + (-p 2 , t) 

+b{-p 1 ,t)[vu} r (p 1 ,p 2 ,t)a{p 2 ,t) + b(-p u f)M r (Pi. p 2 , t)b+{-p 2 , t)}. (9) 

Here and below the vectors pi, p 2 , . . . denote the canonical momenta of fermions 
and ki, k 2 , . . . correspond to the momenta of photons; 

[C»?]^a(Pi>P2;*) = ^a{Pi,t)j' 1 rif3(p 2 ,t)e' r fl , and e 1 ,e 2 are the polarization unit 
vectors, e 3 = k/k and e° = <5°. The photon correlation function is defined as 

F rr ,(k,k',t) = (A+(k,t)A-,(k',t)). (10) 
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With help Eq. ^ we can then obtain the first equation of the BBGKY hierarchy 



F rr i (k, k , t) 



ie(27r) 



-3/2 



]T J d 3 Pl d 3 P2 [ - -4=<^(Pi - P2 - k) 

Q/3 



H^a(Pl' P2> *)KUpi> *)^(P2, } (k, t)) 

M5 Q (P1, P2, t)<0+ (pi, t)b+{-p 2 , (k, *)> 

[w]^a(Pl> P2j t)(&a(-Pl,*)MP2, *)4 _) ( k , *)) 

H^ Q (pi, p 2 , t)<6 (-pi, i)6^(-P2, t)4~ } (k, 0) 
' <5(pi-p 2 + k') 



H^(Pi, P2, t)<4 (Pi, *)a/j(P2, (k', *)) 

+ H^( Pl ,p2,t)(a+(p 1 ,i)6+(-p 2) t)4+ ) (k',i)) 
+ [w]^(Pi,P2, *){6a(-pi,t)a^(p2, t)A { +\k', t)) 

+ [vv} r a (p 1 ,p 2 ,t){b a (-p 1 ,t)b+(-p2,t)A ( r t\k',t)) }.(11) 

3 Truncation procedure 

The simplest decoupling of the hierarchy, 

(a+( Pl ,t)a /3 (p 2 ,i)4 ±) (k,i)) ~ (a+(pi,t)a p (p2,t))(A<&fat)) = 0, (12) 

is not effective. Therefore, we will consider the equation at the second order for 
the annihilation process (the next to last line in Eq. (jlip l and the reverse one 
(the fourth line). We will not write the equations of the second level for these 
correlators completely in view of their awkwardness and will discuss only the 
simplest truncation scheme resulting in correlators of the type 



(a+( Pl , i)a( P2 , t)4±)(k, (k', t)) ~ 

(o+ ( P i , t)a(p 2 , f)) (4±) (k, t)A^ (k', i)) , 



(13) 



which occur in the RPA approximation as well. We ignore other processes 
as, e.g., the vacuum polarization effects contained in correlators of the type 
(a(pi, t)b(— P2, t)Ar \k, t)A^'(k', t)) leading to the polarization functions |5]). 
Ignoring spin effects (see Sect. 2), the approximation p3[) in combination with 
the diagonalization of the photon and fermion correlation functions 



(4+) (k, t)A { r 7 ] (k, t)) = 5 rr ,S(k - k')F r (k, t), 
(P, t)ap(p', t)) = S a0 S(p - p ')/(p, t) 



(14) 
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leads to the following photon KE for zero initial condition 



^(M) = / d3 P It dt'K (p,p -k;t,f)[l + F%tf)] ■ (15) 

[/(P, t') + /(p - k, f ) - 1] cos{/; dr[ W (p, r) + w(p - k, r) - *] }, 

where the nonlocal kernel is K(p,p';t,t') — [vu] r ^ a (p,p';t)[uv] 1 al3 (p / ,p;t'). 
The next step is based on the Markovian approximation that allows to ignore 
the memory effect in the photon distribution, i.e. on r.h.s. of Eq. (|16p we replace 
F(k,t') — > F(k,t). That brings us to the following quadrature formula 

F(k, t) - exp [$(k, t)] - 1 ~ <f(k, t) , (16) 

where the last approximation is valid for the subcritical fields and 

*) = - A ^ d *' d '" / p - k ; *") 

[/(p, O + /(P - k, t") - 1] cos{tf, dr [u(p, r) + w (p - k, r) - fc] }. 

The kernel if(p, p '; i') is a slowly varying function of the momentum argu- 
ments p,p' at fixed t and t'. There are also some complicated fast temporal 
oscillations on this background. For a rough estimate of the effect let us substi- 
tute the kernel K by its average value K (p, p — k; t, t') — > iCo = —5. In addition, 
we neglect the non-Markovian effect in the fermion distribution functions. Since 
the main part of fermions is created from vacuum with small momenta they can 
be neglected due to the momentum and field (E <§C E c ) dependence in the high 
frequency factor on the r.h.s. of Eq. (|17p . We obtain the result 

where 5 — 2m — k is the frequency mismatch and n(t) = 2 J d 3 pf(p, t)/(2ir) 3 is 
the pair density with the factor 2 from the spin degeneracy. Thus, in the optical 
region k -C m the distribution (JT7J) gives F(k) ~ 1/k, which is characteristic 
for the flicker noise (see, e.g., [10]). In the high frequency region k ~ m the 
Markovian approximation is not justified. Here, a more detailed investigation 
is necessary. 



Summary 

Our main result in Eq. (|17p defines the frequency dependence of the photon 
distribution by the factor 1/fc. That this is a multi-photon process can be 
seen when one identifies the frequency mismatch with the energy Nk of the 
photon system which is necessary for the energy conservation in the one-photon 
annihilation process: we obtain N ~ 2m/k. For optical lasers this is a huge 
number and therefore such kind of events are very rare. This conclusion about 
the role of multi-photon processes is conform with the analysis of the absorption 
coefficient of the electron - positron plasma created from the vacuum in the 
infrared region [19] . 
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